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Abstract 

We study the evolution of primordial magnetic fields in an expanding cosmic plasma. For this purpose we present 
a comprehensive theoretical model to consider the evolution of MHD turbulence that can be used over a wide range 
of physical conditions, including cosmological and astrophysical applications. We model different types of decaying 
cosmic MHD turbulence in the expanding universe and characterize the large-scale magnetic fields in such a medium. 
Direct numerical simulations of freely decaying MHD turbulence are performed for different magnetogenesis scenar¬ 
ios: magnetic fields generated during cosmic inflation as well as electroweak and QCD phase transitions in the early 
universe. 

Magnetic fields and fluid motions are strongly coupled due to the high Reynolds number in the early universe. 
Hence, we abandon the simple adiabatic dilution model to estimate magnetic held amplitudes in the expanding uni¬ 
verse and include turbulent mixing effects on the large-scale magnetic held evolution. Numerical simulations have 
been carried out for non-helical and helical magnetic field configurations. The numerical results show the possibil¬ 
ity of inverse transfer of energy in magnetically dominated non-helical MHD turbulence. On the other hand, decay 
properties of helical turbulence depend on whether the turbulent magnetic held is in a weakly or a fully helical state. 

Our results show that primordial magnetic fields can be considered as a seed for the observed large-scale magnetic 
fields in galaxies and clusters. Bounds on the magnetic held strength are obtained and are consistent with the upper 
and lower limits set by observations of extragalactic magnetic fields. 


‘Electronic address: tinatin@andrew.cmu.edu 
^Electronic address: Axel.Brandenburg@Colorado.edu 
^Electronic address: aleko@tevza.org 
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I. INTRODUCTION 


Understanding the origin and evolution of cosmic magnetism is one of the challenging questions of mod¬ 
ern astrophysics. The major questions include theoretical as well as observational aspects of the problem: 
when and how was the cosmic magnetic field generated? How did it evolve during the expansion of the 
universe? What are modern observational constraints on the magnetic fields af large scales? Are magnefic 
fields observed af galactic and exfragalacfic scales of cosmological or asfrophysical origin? The fypes of 
furbulence considered here are characferized by a sfrong random inifial magnefic field. The interacfion wifh 
fhe velocify field leads foo inverse specfral fransfer foward large scales fhaf is unknown in non-magnefic 
furbulence.^ 

The goal is fo identify imporfanf properties of cosmic magnefic furbulence in fhe expanding universe. 
Properties of decaying MHD furbulence in primordial plasma link magnefogenesis scenarios operating in 
fhe early universe wifh fhe consfrainfs on fhe large-scale magnefic fields sef by presenf observations. Hence, 
sfudying fhe magnefic field evolution, we can identify likely magnefogenesis scenarios responsible for ex- 
cifing seed fields in fhe early universe and exclude unlikely ones using consfrainfs sef by modern or fufure 
observations. 

The problem of cosmological magnefogenesis is guided by recenf observafions of large-scale magnefic 
fields. Indeed, galaxies are known fo have magnefic fields fhaf are partly coherenf on fhe scale of fhe galaxy 
wifh field sfrengfhs reaching 10“® Gauss (G) (see Refs. Ilj-S] and references fherein). These magnefic 
fields are fhe resulf of amplificafion of inifial weak seed fields of unknown nafure. Moreover, if is now clear 
fhaf /iG-sfrengfh magnefic fields were already presenf in normal galaxies (like our Milky Way) when fhe 
universe was less fhan half of ifs presenf age ITHa]. This poses sfrong limifs on fhe seed magnefic field 
sfrengfh and ifs amplificafion fimescale. 

From a fheorefical poinf of view fhere are fwo scenarios fhaf can lead fo fhe generation of magnefic 
fields af exfragalacfic scales ilOll : a boffom-up (asfrophysical) scenario, where fhe seed field is fypically 
very weak and fhe observed large-scale magnefic field is fransported from local sources wifhin galaxies fo 
larger scales illH . and a fop-down (cosmological) scenario where a significanf seed field is generafed prior 
fo galaxy formation in fhe early universe on scales fhaf are large af fhe presenf fime il2ll . The major fheme 
of fhis review is fo discuss fhe evolution, sfrucfure, and effecls of cosmic magnefic fields wifh fhe goal fo 
beffer undersfand ifs origin and observational signafures. 

We will briefly discuss cosmic magnefohydrodynamic (MHD) furbulence in order fo undersfand fhe 
magnefic field evolution. MHD furbulence in fhe confexf of asfrophysical plasma processes has been sfudied 
for a long fime. On fhe ofher hand, fhe effecfs of MHD furbulence in cosmological confexfs has received 
affenfion only in recenf years il3n . Simulafions show fhaf fhe kinefic energy of furbulenf motions in galaxy 
clusfers can be as large as 5-10% of fhe fhermal energy densify il4l] . This can influence fhe physics of 
clusfers ilSll . and af leasf should be modeled correcfly when performing large-scale simulations il6l42in . 
Turbulenf mofions can also affecf cosmological phase fransifions; see Refs. i22M24l] and references fherein). 
Turbulence can be generafed by a small inifial cosmological magnefic fields. Undersfanding mechanisms for 
excifing primordial furbulence is an imporfanf goal. We argue fhaf even if fhe fofal energy densify presenf in 
furbulence is small, ifs effecfs mighf be subsfanfial because of fhe sfrongly nonlinear nafure of fhe relevanf 
physical processes. 


' The paper is based on the presentation by Tina Kahniashvili Cosmic Magnetic Fields: Origin, Evolution, and Signatures at the 
Turbulent Mixing and Beyond Workshop 2014 'Mixing in Rapidly Changing Environments - Probing Matter at the Extremes’. 
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Recent important observations Refs. 12514331] (also Ref. 1341] for recent study, and Ref. ll35h for discus¬ 
sions on possible uncertainties in the measurements of blazar spectra), suggest the existence of magnetic 
fields in the universe at scales large enough to suggest a primordial origin ilOl] . This result is robust to po¬ 
tential plasma instabilities of the two-stream family 13614381] . Prior to these observations, there existed only 
upper limits of the order of a few nG for the intergalactic magnetic field. These were obtained through Fara¬ 
day rotation of the cosmic microwave b ackgrou nd (CMB) polarization plane 113914491] and Faraday rotation 
of polarized emission of distant quasars i5f]l453l] . Other tests to derive upper limits on large-scale correlated 
magnetic fields are based on fheir effecf on fhe CM B (se e Ref. 1541] and references fherein), (SJ- 821]. CMB 
disforfions |83-92], fhe broken isofropy l imifs , J^ - 102 ]. big bang nucleosynfhesis (BBN) dafa 1103- 105 ]. 
or large scale sfrucfure (LSS) formalion il06l4l28l] . The lower limit on fhe infergal acfic magnetic field in 
voids of order 10“^® G on 1 Mpc scales is a puzzle of modern asfrophysics (see Ref. il291] ). and could very 
well be fhe resulf of fhe amplificalion of a primordial cosmological field B l|] . 

In whaf follows we review recenf efforls which include fhe pioneering sfudies of primordial magnefic 
field evolufion fhrough cosmological phase fransifions; see Refs. il30l4l40l] . The decay of cosmic magnefic 
field in fhe universe has been analyzed fhrough numerical simulafions of decaying MHD furbulence. Major 
findin gs include (i) fhe possibilify of fhe inverse Iransfer of non-h elical causally generafed magnefic fields 
1 13.51] ; (ii) fasl growfh of vorficify in fhe magnetized universe il32l] ; (iii) growfh of heli cal strucfures a t larg e 
scales for partially helical magnetic fields generated at cosmological phase-transitions il33Lll34lll36Lll37l] . 
and more intere stingly th e absence of the inverse cascade for inflation-generated fully or partially helical 
magnetic fields 1 139 . 140 1. 


II. MODELING MHD TURBULENCE IN THE UNIVERSE 


The orig in of the cosmic magnetic field has been discussed for decades, starting with Enrico Fermi’s pa¬ 
per of 1949 il4in . The approach presented below is novel in several ways, (i) Primordial magnetic fields are 
generally analyzed in the “frozen-in” approximation due to a high conductivity of cosmic plasma, when the 
magnetic field evolves only due to the dilution of field lines as the universe expands. In contrast, we account 
for the actual coupling between the magnetic field and the cosmic plasma, which leads to major differences 
with the frozen-in approximation at some epochs, (ii) Much work on MHD turbulence is focused on specific 
astrophysical objects (such as galaxies, clusters, interstellar medium, or stellar magnetosphere). Instead we 
have developed a comprehensive theoretical framework to consider the evolution of MHD turbulence over 
a wide range of physical conditions, beyond any specific application, (iii) Cosmic MHD turbulence is 
usually studied within one of two limiting cases, the viscous (optically thick) or free-streaming (optically 
thin) regimes. These two regimes differ in the form of viscous or drag forces. Realistic turbulent behavior 
is somewhere in between these two limits, and the numerical simulations have the capability to describe 
adequately a smooth transition between these two regimes. 

As noted above, several astrophysical observations show the presence of a large-scale correlated mag¬ 
netic fields in the universe. The recent study by Dolag et al. Bin concludes that these magnetic fields are 
most likely seeded by a field of primordial origin. In fact, many different mechanisms of cosmological seed 
magnetic field generation have been p roposed. Some of these employ symmetry breaking during phase 
transitions (e.g. electroweak or QCD) il42l4l62l] . On the other hand, if the magnetic field originated dur¬ 
ing a cosmological phase transition, its configuration is strongly limited by causality 11631] : the correlation 
length of the magnetic field cannot exceed the Hubble horizon at the moment of field generation. The 
causality condition combined with the divergence-free field condition implies a magnetic energy spectrum 
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at large sc ales EM{k) oc 11641] (the so-called Batchelor spectrum il65h l. Recent numerical simulations 
1 13ll4l35n confirm that cosmological turbulence produces a Batchelor spectrum completely independently 
of initial conditions present in the cosmic plasma. Combining this causal spectrum with the requirement 
that the total energy density of the magnetic field be less fhan 10% of fhe radiation energy density (fo be 
consisfenf wifh sfandard BBN) leads fo a sfrong limif on fhe smoothed amplitude of the magnetic field 
at large scales of the order of 10“^® to 10“^® Gauss at 1 Mpc il66ll . although the effective value of th e 
magnetic field derived through its total energy density is high enough, of the order of 10“® Gauss lll3fll] . 
(Note that this argument does not account for further evolution of the magnetic field in MHD turbulence). 
Taking into account that the magnetic field effects are mostly determined by effective values (i.e. the total 
energy density), and noticing that the extremely low limits at large scales of causal fields are consequences 
of normalization (smoothing procedure), the upper b ound s have been re-determined in terms of the effec¬ 


tive strength of the magnetic field; see Refs. 11441 .1461 1 12211 . The BBN limits have also been re- analy zed by 


accounting for the MHD evolution of toy magnetic fields throughout expansion of the universe 113411 . 

Ou r particul ar interest lies in helical magnetic fields that can be generated in the early universe; see 
Refs. lll67l4l79t] and references therein. There are two main motivations for considering helical seed mag¬ 
netic fields: (i) the presence of h elical magnetic fields in the early universe can be related to the lepto- 
and baryogenesis problems 1 II 8 OI 1 ; (ii) it she ds light on the evolution of helical magnetic fields in stellar 
magnetospheres, AGNs, and voids il8ll.ll82ll . 

An exception to the Batchelor spectrum (spectral index n = 4) is the possibility of inflationary magne¬ 
togenesis, in which the spectral index of th e magnetic field could be less than -|-1, and the simplest option 
is a scale-invariant spectrum with re —)• — 1 11831420711 . Inflation-generated magneti c fie l d sce narios should 
be considered with some caution due to the possibility of significant bac krea ction 11208142 llh . which is not 
an issue for the phenomenological, effective classical model; see Ref. 121211 and references therein. The 
first simulations describing the inflation-generated magnetic field coupled to the primordial plasma sug- 
geste d that the presence of an initial magnetic field leads to large-scale turbulent motions in the rest plasma 
1 13211 . Ongoing research consists in the study of inflation-generated he lical mag netic field (with a scale- 
invariant k~^ spectrum) evolution during the expansion of the universe il38l4l40ll . Simulations show that 
inflation-generated magnetic fields retain information about initial conditions. In other words, they decay 
very slowly when compared with phase transition-generated fields. Magnetic fields are almost “frozen-in” 
the primordial plasma at large scales, where causality allows interaction only at scales smaller than the 
Hubble horizon and they correspondingly retain their initial spectral shape. On the other hand, within the 
causal horizon, the magnetic seed field interacts with cosmic plasma leading to the excitation of kinetic 
motions (turbulent velocities); see below. We also discuss an alterna tive appro ach where cosmic magnetic 
fields originate during the late stages of the evolution of the universe BlOfil. 121311 . In this case the correlation 
length is strongly limited by the causality requirement. Due to the sharp spectral shape at large scales, the 
magnetic field amplitude might be low. The simplest astrophysical magnetogen esis mech anism invokes the 
ejection of magnetic flux from compact systems such as AGNs or supernovae 12141.121511 . In this scenario 
the generation of a strong magnetic field is ensured by its extremely fast generation due to rapid rotation 
of the object lllll . Other mechanisms are based on the generation of a small see d by plasma processes 
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III. THE MODEL 


As mentioned above, we focus on the magnetic field evolution during the expansion of the universe 
from the moment of magnetic field generafion until today. Over fhis lengfhy period, fhe magnetic field is 
affected by difterenf physical processes fhaf resulf in amplificalion as well as damping: fhe complexifies of 
fhe problem are due fo fhe sfrong coupling befween magnetic field and furbulenf mofions. Firsf, fo accounf 
for fh e cosmological expansion we musf reformulate the MHD equations in terms of comoving quantities 
1 2201] ■ Specific epochs mosf relevanf fo fhe final configuralion of fhe primordial magnetic field are related 
to cosm olog ical phase transitions, neutr i no decoup ling, nucleosynthesis, recombination, and reionization; 
see lie, 221] for reviews and Refs. 1 110 . 222 - 227 1. 

In the following, we discuss numerical simulations performed with t he P ENCIL CODE. This public 
MHD code (https : //github . com/pencil-code) (see also Ref. 122811 1 is particularly well suited 
for simulating turbulence owing to its high spatial (sixth order) and temporal (third order) accuracy, while 
still taking advantage of the finite difference in terms of speed and straightforward parallelization. Recent 
results fro m the PEN CIL CODE include MHD turbulence simulations at the electroweak or QCD phase 
transitions 1 131 - 135 1. 

a. Numerical Technique. By default, the PENCIL CODE solves the MHD equations for the logarith¬ 
mic density In p, flow velocify v, and fhe magnefic vector pofenfial A as follows: 


D , 

ft;'"" 


D 

Drj 

d 


-V = 


dr] 


A = 


—V • V , 

J X B - c^Vlnp-b /vise , 
V X B f]\4 -|- AV^A . 


( 1 ) 

( 2 ) 

(3) 


Here, p is fhe conformal lime and D/Dr/ = O/c)?/ -|- rr • V is fhe adveclive lime derivative, /vise = 
zr -|-2 VV • is fhe compressible viscous force for conslanf v, Gi = 2 SjjVjln/ 9 , and 

~ “ \^ijVk,k is the traceless rate-of-strain tensor. The pressure is given by p = pel, 

where = 1/\/3 is the speed of sound in the case of an ultra-relativistic gas, and J = V x S/dvr is the 
current density. 

In our simulations we use a vanishing magnetic forcing term /m = 0 everywhere, except for the purpose 
of producing initial conditions, as explained below. 

b. Initial Conditions. To produce initial conditions, we run the simulation for a short time (Af 
O.SAi / Cs) with a random (in time) (5-correlated magnetic force /m in Eq. (3). The forcing term is composed 
of plane monochromatic waves pointing randomly in all possible directions with an average wavenumber 
feo and fractional helicity (/m • V x fu) / (^o/m) =2a/{l + tj^). Here a is the parameter characterizing 
the initial forcing. Initial conditions for the magnetic and velocity fields produced from such a procedure 
have fhe advantage of being turbulent, still self-consistent solutions of the MHD equations. 

c. Ejfective Magnetic Field Characteristics. A magnetic field generaf ed during pha se fransifions 

fhrough any magnefogenesis scenario should satisfy fhe causalify condition lll63L Il64l. 122011 . The maxi¬ 
mal correlation lengfh ^max of a causally generated primordial magnetic field should be shorter fhen fhe 
Hubble radius af fhe time of generafion, We define fhe paramefer 7 = Cmax/2?/^ < 1, which can 

describe fhe number of primordial magnefic field bubbles inside fhe Hubble radius, and fhus N oc 7 “^. 
To accounf for fhe universe expansion we use comoving lengfh, which is measured today and corresponds 
fo fhe Hubble radius af fhe momenf of magnefic field generafion. Comoving lengfh should be inversely 
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proportional to the phase transition temperature (T*): 


Ah* = 5.8 X 10"^° Mpc 


I'lOOGeV^ 




(4) 


For the QCD phase transition ( 5 * = 15 and T* = 0.15 GeV), the comoving length equals 0.5 pc, while 
for the electroweak phase transition (p* = 100 and T* = 100 GeV) it should be equal to 6 x 10“^ pc. In all 
cases, the correlation length of the primordial magnetic field should not exceed the comoving value of the 
Hubble radius: ^max < Ah- Obviously, the latter condition accounts only for the increase of the correlation 
length due the expansion of the universe, and does not account for the effects of cosmic MHD turbulence 
(free decay or an inverse cascade in the case if primordial magnetic fields have nonzero helicity). Note, that 
the number of bubbles inside the Hubble radius is around 6 (7 ~ 0.15) for the QCD phase transition and 
around 100 (7 ~ 0.01) for the electroweak phase transition. Thus, the maximal correlation length for the 
QCD and electroweak phase transitions should be 0.08 pc and 6 x 10“® pc, respectively. On the other hand, 
the correlation length is unlimited in the case of inflation-generated magnetic fields. 

As mentioned above, the primordial magnetic field contributes to the relativistic component and thus the 
total energy density of the primordial magnetic field where on is the scale factor during nucleosyn¬ 

thesis, is limited by the BBN bound: it cannot exceed 10% of the radiation energy density /9rad(<^N)- It is 
straightforward to see that the maximal value of the effective magnetic field defined thr'ough the total mag¬ 
netic energy does not depends on the temperature at the moment of generation (T*), and depends weakly on 
the relativistic degrees of freedom (7) at the moment of the magnetic field generation. 

The dominant contribution to the magnetic field energy density comes from the given length-scale, 
the so-called integral scale, where the magnetic field strength reaches its maximum. Thus, when dealing 
with phase transition-generated magnetic fields, we adopt the following idealizing approximation: we 
generate initial conditions for freely decaying turbulence simulations by running a numerical simulation of 
forced MHD equations for a short time interval. The external electromagnetic force, intended to generate 
a turbulent state, is introduced in the form of 6 functions that peak at a characteristic wavenumber, kg = 
27r/^Q^. This yields random magnetic fields with correlation length ^o- Thus, the magnetic field strength 
at the characteristic length scale is = s/Sttpb- The characteristic length scale of the initial turbulent 

state ^0 is set by the size of the largest magnetic eddies, because the primordial magnetic field evolves with 
the primordial MHD turbulence. In this approach, the characteristic length scale of the magnetic field is set 
by the bubble size of the phase transition when magnetogenesis occurs. 

d. Magnetic Field Spectrum. The interaction between magnetic field and plasma gives rise to kinetic 
motions, and the turbulent backreaction results in spreading of the spectral energy density of magnetic field 
over a range of wavenumbers. At scales longer then the integral scale of the turbulence (small wavenum¬ 
bers), the spectral energy density develops into the form of a power law Em = Afc”, where A is a normaliza¬ 
tion constant, and n is the spectral index. The spectrum of the turbulent magnetic field can be determined by 
the spectral expansion of the two-point correlation function of the magnetic field (i?j(x)i?j (x -|- r)), whose 
Fourier transform with respect to r gives the spectral function 

-Gj (k) - (k) • G) 


Here, Pij (k) = 5ij—kikj/Siji is the antisymmetric tensor, and H^{k) is the magn etic helicity spectrum. 
In this case, a white noise spectrum corresponds to the spectral index n = 2 lll63n . while the Batchelor 
spectrum corresponds to the spectral index n = 4 lll64l] . The power law of large-scale MHD turbulence 
spectrum extends down to the integral scale ^m, which is itself a time-dependent quantity throughout the 


7 








turbulence decay process. At scales shorter than the integral scale the spectral energy density of the magnetic 
field decreases rapidly due to the combined action of turbulent decay and viscous damping. 

e. Decay Laws. The correlation length of the turbulent magnetic field evolves in lime during Ihe free 
decay of lurbulence. We may describe Ihe decay laws of Ihe magnetic correlalion lenglh ^Miv) Ihe 
speclral energy density SMif]) using Iwo power law indices and ue- 

= ^Mivo) f , (6) 

\VoJ 


SMiv) = 


riE 


(V) 


The speclral energy densify of Ihe primordial MHD lurbulence speclrum can be splil info ils large-scale and 
shorl-scale componenls, above and below Ihe time dependenl integral scale: 


EM{k,r]) = Eo{rj) 


when k < kj{r]) 
k~^/^ when k > ki{r]) 


( 8 ) 


where k = k/kj and ki{r]) = 27r/^M(^)- Hence, Eqs. (6) and (7) can be used to describe the time evolution 
of spectral amplitude of magnetic field Eq (ry) for a given turbulent spectrum: 


Eo{ri) 




(9) 


IV. RESULTS 


A. The Inverse Transfer for Non-Helical Fields 


The inverse cascade is by now a well-known effect in helical magnetic turbulence 12291] . One of the 
remarkable results is the presence of non-helical inverse transfer for magnetically dominated (causally 
generated) MHD turbulence; see Fig. [T] where we show spectral energy transfer rates, which demonstrate 
that the inverse transfer is about half as strong as with helicity. However, in both cases the magn etic g ain at 
large scales results from velocity at similar scales interacting with sma ller-scale magnetic fields 113511 . This 
resull has nol been emphasized in previous sludie s, se e Ref s. i230Ll23in and references Iherein, and has now 
been confirmed by independenf research groups 12321423411 . 

Recenf high resolution simulafions wifh differenl magnefic Prandfl numbers PrM = i^/A 113511 have 
shown a clear k~‘^ speclrum in Ihe inerfial range. This is Ihe firsl example of fully isolropic magnelically 
dominaled MHD lurbulence (governed by fh e phase Iransilion-generaled magnefic fields) exhibiting whal 
we argue lo be weak lurbulence scaling 123511 . On Ihe ofher hand, Ihe Kolmo gorov scaling k~^^^ has been 
recovered for Ihe case of Ihe inflalion-generaled helical magnefic fields 113911 . 


B. Inflation Generated Magnetic Fields 

Magnetic fields generated during the inflation should be affected by cosmological phase transitions 
occurring at later times during the expansion of the universe. In this case, a separate study of the imprint of 
phase transitions on cosmic magnetic fields is needed. For this purpose we adopt a general approach that can 














FIG. 1: Spectral transfer function T^pq, (a) as a function of k and summed over all p and q, (b) as a function of p and 
q for /c//ci = 4, and (c) as a function of k and q for p/ki = 4. The dashed line in (a) and the insets in (b) and ( c) sh ow 
the corresponding case for a direct numerical simulations with helicity; both for PrM = 1- See Fig. 3 of Ref. 1135 1. 


be applied to both, QCD and electroweak phase transitions. In each case, turbulence forcing is determined 
by the phase transition bubble size. Rapid phase transitions generate turbulence, which then decays slowly 
at large scales. In contrast to previous studies, the inflation-generated magnetic field is not frozen into the 
cosmic pla sma. Turb ulence is generated during a short forcing period, which then is followed by slow decay 
(see Refs. il3ll. 113211 for details). Recent simulations showed an increasing characteristic length scale of 
the velocity field and the establishment of a k'^ (white noise) spectrum at large sc ales. This increase of 
vorticity perturbations occurs until it reaches equipartition with the magnetic field 113411 . Figure |2] shows 
the evolution of kinetic and magnetic field spectra from those simulations. Numerical results show that 
inflation-generated magnetic fields are not significantly modified at large scales by their coupling to the 
plasma during a cosmological phase transition. The coupling of cosmic magnetic field with the phase 
transition-generated fluid turbulence leads to deviations of the magnetic field spectrum from the initial 
scale-invariant shape only at intermediate scales. 

Ongoing research consists in pursuing high resolution numerical simulations of helical inflation¬ 
generated magnetic field evolution. Such a field, being subject to inflationary expansion, is characterized by 
a scale-invariant spectrum n —)• — 1, and its correlation length can be as large as Hubble horizon today or 
even larger (i.e., even when the total energy density Em is finite, the correlation length cx f dkEM (k)/k 
divergences fo r k —)• 0). In contrast to well known helical magnetic field decay laws il3Lll65Ll226l. |236 - 


238Ll240l4244l] . an absence of the inverse cascade has been found for inflation-generated magnetic fields. 


Furthermore, an unusually slow growth of the correlation length and conservation of helicity has been re¬ 
covered even for the case of partially helical magnetic fields. These unexpected and unknown features 
of magnetic helicity are the result of a substantial turbulent power at large scales and the impossibility of 
the redistribution of helical fields at small wavenumbers (only the forward cascade is possible). A more 
thorough investigation of this phenomenon will be performed through varying initial conditions and basic 
parameters of primordial plasma. 


C. Growth of Helical Structures 


It is long known that the magnetic helicity plays a crucial role in determining the evolution pattern of 
MHD turbulence. Distinct evolution characteristics are known for helical and non-helical fields. In recent 


simulations, a partially helical initial magnetic field was used il33h . assuming a tiny initial magnetic helicity 
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FIG. 2: Magnetic (solid lines) and kinetic (dashed lines) energy spectra in regular time intervals. Re r; 170. The 
magnetic and kinetic spectra at the last time are additionally marked in red and blue, respectively. 



FIG. 3: Visualizations of (upper row) and Vx (lower row) at three times during th e mag netic decay of a weakly 
helical field with cr = 0.03 generated during QCD phase transitions. See Fig. 2 of Ref. lll33ll . 


during the QCD phase transition. It was shown that at late times the resulting field attains the maximally 
allowed magnetic helicity. This result is important since helicity crucially affects the MHD dynamics , 
and has very interesting consequences in astrophysical objects (e.g. galactic magnetic fields i245t4247ll . 
for example). The resulting magnetic field has an amplitude of around 0.04 nG and a correlation length 
of order 20kpc, which (assuming realistic scenarios of amplification 113IH ) serves as a seed for galactic 
magnetic fields. At this point the electroweak phase transition-generated magn etic fields are less promising 
due to a smaller initial corr elation length, but are not completely excluded i248n . in particular for fully 
helical magnetic fields 113411 . 

Magnetic helicity is a crucial factor that affects the evolution of primordial magnetic fields. The evo- 
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FIG. 4: Comparison of (upper row) and In p (lower row) for an inflation-generated magnetic field with a = 1 
(left) and tr = 0.03 (right). 


lution of the primordial magnetic field that has been produced with weak initial magnetic helicity that 
undergoes two consecutive stages. During the first stage, the evolution of a partially helical magnetic field 
spectrum is much similar to that of non-helical magnetic fields, and is sometimes described as a direct 
cascade. At this stage the spectral energy density cascades from large to small scales, where it undergoes 
de-correlation and viscous damping. Magnetic helicity is conserved and hence its fractional value increases 
during the turbulent decay process. The second stage in the primordial magnetic field evolution sets in when 
the turbulent state with maximal helicity is reached. The maximal value of the magnetic helicity that can 
be reached is limited by the realizability condition. Indeed, the conservation of magnetic helicity leads to a 
decay of the magnetic energy density inversely proportional to the correlation length of the turbulence: 

= \nM{'n)\/2£M{ri), ( 10 ) 


where is the minimal correlation length of the turbulent state. He nce, an inverse cascade develops 

during the second stage of the primordial magnetic field evolution. In Ref. il33n . we have studied the 
and for 0.1, and 0.03 in the case of the QCD phase transition. It seems that, especially at 

lower cr, the increase of is slow (~ while ^M{r]) S> However, since magnetic helicity 

conservation implies that £m decreases as ri~^, the minimal correlation length soon reaches 

see Fig. (I5]l. When the correlation length of the turbulent magnetic field reaches a minimal correlation 
length, the turbulence reaches its fully helical state. Then the turbulence decays according to the helical 
turbulent decay laws: ~ and £m ~ Hence, we identify two distinct phases in the MHD 

turbulence evolution: the phase of weakly helical turbulence decay with = 1/2 and ue = —1, and 
the phase oi fully helical turbulent decay with = 2/3 and he = —‘Ij'i). The fully helical case is 
characterized by an inverse cascade whe re Enin) oc CM('n)£M('n ) = const (see Eq. dO]). These results are 
in full agreement with earlier works j 1 3l . 1 1 65112361423 81 124(1124211 . The effective coupling of the primordial 
magnetic fields and cosmic plasma ends at fhe time of recombination. At later stages, primordial magnetic 
fields exhibit much slower developments 122011. 
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FIG. 5: Evolution of turbulent correlation length (solid) and mi nima l correlation length (dashed) of a 

helical state for <7=1 (black), 0.1 (blue) and 0.03 (red). (Fig. 3 of Ref. 1 133 11. 


Knowing the initial values of the turbulent magnetic correlation length the minimal correlation 

length set by the realizability condition, we can calculate the time interval r/fuUy needed for the tur¬ 

bulence to reach its fully helical state during the decay process. Since these two scales approach each other 
as r/t/2^ the result is r/fuUy = Hence, the time interval needed for the development 

of a fully helical state can be calculated using the initial values of turbulent energy and helicity: 

^fully = /'^M- (^^) 

Note that this time increases inversely proportional to the square of the initial helicity of the turbulent 
state T-Lm- Assuming that the initial magnetic helicity is times lower then the maximal helicity in 

the case of strong CP violation during phase transition, we can calculate the time needed for the cosmic 
turbulence to reach its maximally helical state: ?7fuiiy = %/7^- 


V. DISCUSSION 


Let us now discuss our results in the broader context of the workshop ‘Mixing in Rapidly Changing 
Environments—Probing Matter at the Extremes’. Our work has demonstrated a rather generic trend of de¬ 
caying MHD turbulence to displ ay an increase of energy at large length scales. This process is well-k now n 
in helical MHD turbulence 122911 . but to a lesser extent it also occurs in non-helical MHD turbulence 113511 . 
Moreover, if there is small fractional helicity initially, this fraction will increase with time proportional to 
the square root of time. At a certain time, the fractional helicity will be 100%, after which the decay of 
energy slows down and the increase of the correlation length speeds up. 

An increase of energy at large length scales is not a common phenomenon in hydrodynamic turbulence 
and has never been found for passive scalars. This special behavior in MHD is likely to lead to uncon¬ 
ventional mixing properties, although this has not yet been well quantified for decaying MHD turbulence. 
However, for statistically stationary turbulence, an increase of energy at small wavenumbers is usually de- 
scrib ed as non-diffusive turbulent transport, which is particularly well known in mean-field dynamo theory 
1 25(11 . In a sense, this is more reminiscent of what might look like “anti-mixing”. This is to some extent due 
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to the fact that vector fields behave differently from scalar fields. This can in part be due fo fhe presence of 
addifional conservation laws. In parficular, magnefic helicify provides an exfremely powerful consfrainf. 

The significance of sfudying decaying MHD furbulence is manifold. On fhe one hand, our resulfs will 
help fo beffer undersfand fhe nafure of cosmic magnefism and will gain insighf info fhe decay laws of cosmic 
MHD furbulence. On fhe ofher hand, our analysis can be applied nof only fo cosmological scales, buf fo 
molecular clouds or even profoplanefary disks where decaying magnefic furbulence can crucially affecf 
fhe global sfafe or fhe formation of local sfrucfures. The resulfs of our research have fherefore imporfanf 
implicafions in many areas, including fluid dynamics, early universe physics, high energy asfrophysics, 
MHD modeling, and large-scale strucfure formalion in fhe universe. 


VI. CONCLUSIONS 


We have discussed fhe evolution of primordial magnefic fields during fhe expansion of fhe universe 
and have addressed some of fhe observafional signafures. The coupling befween fhe magnefic fields and 
cosmic furbulence leads fo novel resulfs as compared fo previously adopfed frozen-in approximafions when 
magnefic field evolufion was considered solely due fo fhe field line dilufion in fhe expanding universe. For 
fhis purpose, we have developed a comprehensive fheorefical model fo consider fhe evolufion of MHD 
furbulence over a wide range of physical condifions, beyond any specific asfrophysical applicafion. 

Numerical simulations have shown novel effecls in fhe evolufion of magnefic fields in cosmic furbulence. 
If seems fhaf inverse fransfer fhaf normally occurs in helical MHD furbulence, can also fake place for non¬ 
helical magnefic fields if fhe MHD furbulence is magnefically dominafed. In fhis case, large-scale velocify 
perfurbafions power up fhe magnefic field, leading fo subsfanfial increase of magnefic power specfra af large 
scales and corresponding inverse fransfer. 

An analysis of fhe inflation-generated magnetic field is carried ouf fo find ouf how fhey are affecfed 
by fhe cosmic phase fransifions (QCD and elecfroweak phase fransifion). Numerical resulfs show fhaf 
large-scale magnefic fields survive phase fransifions, and fhus, phase fransifions cannof rule ouf inflafionary 
magnefogenesis as fhe source of seed magnefic field in fhe universe. 

On fhe ofher hand, magnefic fields generafed during phase fransifions can have finy helicify. We have 
shown fhaf during fhe evolufion in MHD furbulence magnefic field helicify grows until if reaches maximal 
helical sfafe. If seems fhaf helicify growfh rafe is fasf enough fo reach maximal helicify well before fhe 
epoch of recombinafion, when fhe primordial magnefic field decouples from cosmic furbulence. 

We have used resulfs of high resolufion simulafions fo sef limifs on fhe large-scale magnefic field gener¬ 
afed during early sfages of fhe evolufion of fhe universe (inflafion or cosmological phase fransifions). Indeed, 
if seems fhaf magnefic fields produced during fhis epochs and subsequenfly modified by cosmic MHD furbu¬ 
lence can reach amplifudes similar fo fhe lower bounds of fhe observafional magnetic fields even accounting 
for fhe effecfs of large-scale furbulenf decay as w ell as additional Alfven wave damping. The exfremely low 
values derived for smoofhed magnefic field 116611 do nof imply fhaf fhe effective magnefic field is also small 
in fhe 1 pc-1 kpc range and cannof lead fo fhe observafional signafures in blazar emission specfra. Using 
fhe effective magnefic field approach we obfain resulfs fhaf do nof depend on fhe specific specfral shape 
of fhe magnetized furbulence. Observafional signafures of fhe magnefogenesis during elecfroweak phase 
fransifion can come from fufure observations, if weak magnefic fields wifh G amplifude and 

few pc correlation lengfh are defecfed. While a somewhaf sfronger magnefic field wifh a correlation lengfh 
of fhe order of kpc mighf indicafe fhe presence of QCD phase fransifion magnefogenesis. In furn, magnefic 
fields wifh exfremely large correlation lengfh (IMpc or higher) will indicafe inflafionary magnefogenesis. 
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iog(fM [MPC]) 


FIG. 6: Evolution of the effective magnetic field and correlation length in the case of magnetogenesis at 
electroweak phase transition (green) and QCD phase transition (orange). Arrows indicate the evolutionary path of the 
strength and integral scale of helical and non-helical turbulent magnetic fields during the radiation-dominated era up 
to their final values. Thick solid line(s) show possible pres ent d ay field strengths and integral scales of the magnetic 
field generated during phase transitions (see Fig. 7 of Ref. 111341] '). 
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